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Abstract: We investigate the general properties of Unified Dark Matter (UDM) fluid 
models where the pressure and the energy density are linked by a barotropic equation of 
state (EoS) p = p(p) and the perturbations are adiabatic. The EoS is assumed to admit a 
future attractor that acts as an effective cosmological constant, while asymptotically in the 
past the pressure is negligible. UDM models of the dark sector are appealing because they 
evade the so-called "coincidence problem" and "predict" what can be interpreted as w^e ~ 
— 1, but in general suffer the effects of a non-negligible Jeans scale that wreak havoc in the 
evolution of perturbations, causing a large Integrated Sachs-Wolfe effect and/or changing 
structure formation at small scales. Typically, observational constraints are violated, unless 
the parameters of the UDM model are tuned to make it indistinguishable from ACDM. 
Here we show how this problem can be avoided, studying in detail the functional form of 
the Jeans scale in adiabatic UDM perturbations and introducing a class of models with a 
fast transition between an early Einstein-de Sitter CDM-like era and a later ACDM-like 
phase. If the transition is fast enough, these models may exhibit satisfactory structure 
formation and CMB fluctuations. To consider a concrete case, we introduce a toy UDM 
model and show that it can predict CMB and matter power spectra that are in agreement 
with observations for a wide range of parameter values. 
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1. Introduction 

In the last three decades the flat ACDM model |fl, [| has emerged as the standard "con- 
cordance" ||, |j] model of cosmology. It assumes General Relativity (GR) as the correct 
theory of gravity, with two unknown components dominating the dynamics of the late 
Universe: i ) a cold collisionless Cold Dark Matter (CDM) describing some weakly inter- 
acting particles, responsible for structure formation, ii ) a cosmological constant A p, |6| 
making up the balance to make the Universe spatially flat and driving the observed cosmic 
acceleration M, ||, P, p!0| . The main alternative to the cosmological constant is a more 



general dynamic component called Dark Energy (DE) |Tl|, 12, [Oj. Many independent ob- 
servations support both the existence of a CDM component and that of a separate DE 

@, m, u, ii, n m n n n. 



Early proposals [|, § of the ACDM model were adding A to CDM in an attempt to 
conciliate in the simplest possible way the emerging inflationary paradigm, which requires 
a spatially flat Universe and an almost scale-invariant spectrum of perturbations, with the 
observed low density of matter. It should however be recognised that, while some form of 
CDM is independently expected to exist within any modification of the Standard Model 
of high energy physics, the really compelling reason to postulate the existence of DE has 
been the cosmic acceleration measured in the last decade |7], ||, ^, [[(], |17|, 18|, 1£, ^0], 21]. It 



is mainly for this reason that it is worth investigating the hypothesis that CDM and DE 
are the two faces of a single Unified Dark Matter (UDM) component, thereby also avoiding 



the so-called "coincidence problem" 23] 



Other attempts to explain the observed acceleration also exist, most notably by as- 
suming a gravity theory other than GR, or an interaction between DM and DE (see e.g. 
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in], H H. H 13 and HI Hi' IH HI H H> II)- In this p a P er however we focus ° n 

UDM models, where this single matter component provides an explanation for structure 
formation and cosmic acceleration. 

In general, in the ACDM model or in most models with DM and DE, the CDM 
component is free to form structures at all scales, with DE only affecting the general 
overall expansion [11, [T^, 13 1. Instead, a general feature of UDM models is the possible 
appearance of an effective sound speed, which may become significantly different from 
zero during the Universe evolution, then corresponding in general to the appearance of a 
Jeans length (i.e. a sound horizon) below which the dark fluid does not cluster (e.g. see 
i. Moreover, the presence of a non-negligible speed of sound can modify the 



evolution of the gravitational potential, producing a strong Integrated Sachs Wolfe (ISW) 
effect [|3"6f| . Therefore, in UDM models it is crucial to study the evolution of the effective 
speed of sound and that of the Jeans length. 

Several adiabatic fluid models and models based on non canonical kinetic Lagrangians 



have been investigated in the literature. For example, the generalised Chaplygin gas [38 
39|, HI (see also gl]> g |£| g|, g|, g7], £§ g g ||), the Scherrer §§ and generalised 
3|, the single dark perfect fluid with "affine" 2-parameter barotropic 

~|]) and the 



Scherrer solutions 



equation of state (see [54, 37] and the corresponding scalar field models 
homogeneous scalar field deduced from the galactic halo space-time [56, 57]. In general, in 
order for UDM models to have a background evolution that fits observations and a very 
small speed of sound, a severe fine-tuning of their parameters is necessary (see for example 
|37|, (|7|, H, H, H> HI' HI)- Fm ally, one could also easily reinterpret UDM models 
based on a scalar field Lagrangian in terms of generally non-adiabatic fluids ^] (see 
also |53| , |62|). For these models the effective speed of sound, which remains defined in the 
context of linear perturbation theory, is not the same as the adiabatic speed of sound (see 
pq] , [pjf and |64|]). In |52] a reconstruction technique is devised for the Lagrangian, which 
allows to find models where the effective speed of sound is small enough, such that the 
A;-essence scalar field can cluster (see also ]65|]). 

In the present paper we investigate the possibility of constructing adiabatic UDM 
models where the Jeans length is very small, even when the speed of sound is not negligible. 
In particular, our study is focused on models that admit an effective cosmological constant 
and that are characterised by a short period during which the effective speed of sound varies 
significantly from zero. This allows a fast transition between an early matter dominated 
era, which is indistinguishable from an Einstein-de Sitter model, and a more recent epoch 
whose dynamics, background and perturbative, are very close to that of a standard ACDM 
model. 

To consider a concrete example, we introduce a 3-parameter class of toy UDM adiabatic 
models with fast transition. One of the parameters is the effective cosmological constant 
p\ or, equivalently, the corresponding density parameter Oa; the other two are p s and 
Pt, respectively regulating how fast the transition is and the redshift of the transition. 
Studying the Jeans scale in these models we find an approximate analytical relation that 
sets a constraint on these two parameters, a sufficient condition that p s and pt have to 
satisfy in order for the models to be minimally viable. This relation can be used to fix p s 
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for any given p t : in this case, with respect to a fiat ACDM, in practice our models have one 
single extra parameter. With the help of this relation we establish our main result: if the 
fast transition takes place early enough, at a redshift z > 2 when the effective cosmological 
constant is still subdominant, then the predicted background evolution, Cosmic Microwave 
Background (CMB) anisotropy and linear matter power spectrum are in agreement with 
observations for a broad range of parameter values. In practice, in our toy models the 
predicted CMB and matter power spectra do not display significant differences from those 
computed in the ACDM model, because the Jeans length remains small at all times, except 
for negligibly short periods, even if during the fast transition the speed of sound can be 
large. In other words, this kind of adiabatic UDM models evade the "no-go theorem" of 
Sandvik et al [50] who, studying the generalised Chaplygin gas UDM models, showed that 
this broad class must have an almost constant negative pressure at all times in order to 
satisfy observational constraints, making these models in practice indistinguishable from 
the ACDM model (see also @)- 

The paper is organised as follows: in section ^ we introduce the basic equations describ- 
ing the background and the perturbative evolution. In section ^ we use the pressure-density 
plane to analyse the properties that a general barotropic UDM model has to fulfil in order 
to be viable. In section |] we introduce our toy UDM model with fast transition and study 
its background evolution, comparing it to a ACDM. In section [| we analyse the properties 
of perturbations in this model, focusing on the the evolution of the effective speed of sound 
and that of the Jeans length during the transition. Then, using the CAMB code ||6q |, in 
section |6] we compute the CMB and the matter power spectra. Finally, section is devoted 
to our conclusions. 



2. Background and perturbative equations for a UDM model 

We assume a spatially flat Friedmann-Lemaitre-Robertson- Walker (FLRW) cosmology. 
The metric then is ds 2 = —dt 2 + a 2 (t)&ijdx l dx :) , where t is the cosmic time, a(t) is the 
scale factor and 5ij is the Kronecker delta. The total stress-energy tensor is that of a per- 
fect fluid: T^y = (p + p) u^Uy + pg^ u , where p and p are, respectively, the energy density 
and the pressure of the fluid, while is its four-velocity. Starting from these assumptions, 
and choosing units such that 8ttG = c = 1, Einstein equations imply the Friedmann and 
Raychaudhuri equations: 

- = ~(p + 3p), (2.2) 
a b 

where H = a/a is the Hubble expansion scalar and the dot denotes derivative with respect 
to the cosmic time. Assuming that the energy density of the radiation is negligible at the 
times of interest, and disregarding also the small baryonic component, p and p represent 
the energy density and the pressure of the UDM component. 
The energy conservation equation is: 

p = -3H (p + p) = -3Hp (1 + w) , (2.3) 
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where w := p/p is the equation of state (hereafter EoS) "parameter". In this paper we 
investigate a class of UDM models based on a barotropic EoS p = p(p), i.e. those models 
for which the pressure is function of the density only (see e.g. |37j] and |68| for a discussion 
of general properties of barotropic fluids as dark components). In this case, if the EoS 
allows the value w = — 1, the barotropic fluid admits an effective cosmological constant 
energy density, i.e. a fixed point of Eq. ( |2.3| ) (67], [6^, [54|], which we will denote as p\. Under 
very reasonable conditions (see the discussion below) this effective cosmological constant 
is unavoidable for barotropic fluids 1 |37|, |SSf] - 

In order to properly describe the dynamics of the fluid we must consider the back- 
ground EoS as well as the speed of sound which regulates the growth of fluid perturbations 
on different cosmological scales. In the following we shall confine our study to the simplest 
hypothesis that the EoS remains of the barotropic form p = p(p) when we allow for pertur- 
bations: in this case our models will be adiabatic, and the effective and adiabatic speeds 
of sound will coincide, see e.g. 170, 71, 72, E^. Other choices for the perturbed spacetime 



are possible, see [37] for a recent practical example. 

Let us consider small perturbations of the FLRW metric in the longitudinal gauge, 
using conformal time r\: ds 2 = —a 2 {rf) [(1 + 2$) drj 2 — (1 — 2$) 5ijdx l dx^~\ , where $ is the 
gravitational potential. 

Defining 



2$ 



u := 



(2.4) 



Vp + p 

and linearising the 0-0 and 0-i components of Einstein equations, for a plane-wave per- 
turbation u oc exp (ik • x) one obtains the following second order differential equation 



|7J, |8 



u" + k 2 c 2 u - —u = , 



(2.5) 



where the prime is the derivative with respect to the conformal time rj, cj is the effective 
speed of sound and 

;(! + *), 



P 

3(p + p) 

with z the redshift, 1 + z = a -1 . In general, the adiabatic speed of sound is c 2 d 
for an adiabatic fluid c 2 = c^ d . 

Starting from Eq. (|2.5|), let us define the squared Jeans wave number [36 



(2.6) 

p'/p'; 



c 2 e 



Its reciprocal defines the squared Jeans length: Aj := a 2 /k 2 . 



(2.7) 



1 Since this effective cosmological constant trivially satisfies an energy conservation equation ( |2.3| ) on its 
own, a fluid admitting an effective pa is always equivalent to two separate components, namely pa itself 
and an "aether" fluid, see |67]] and ps} |. Obviously, this is more general; for instance, a scalar field with a 
potential admitting a non vanishing minimum Vo = V((j>o) 7^ at - say - cj>o, is equivalent to a cosmological 
constant p\ = Vo and a scalar field in a potential V = V — Vo- 
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There are two regimes of evolution. If k 2 3> k 2 and the speed of sound is slowly 
varying, then the solution of Eq. ( |2.5[ ) is 



u ~ — — exp ( ±ik I c s drj \ , (2.£ 



where C is an appropriate integration constant 2 . On these scales, smaller than the Jeans 
length, the gravitational potential oscillates and decays in time, with observable effects on 
both the CMB and the matter power spectra |]36| . 

9"/0, with general solution 

drj 



For large scale perturbations, when k 2 <C k 2 , Eq. (|2.5| ) can be rewritten as u" /u 



u^ Kl e + K 2 9 J . (2.9) 

In this large scale limit the evolution of the gravitational potential <1? depends only on the 
background evolution, encoded in 8, i.e. it is the same for all k modes. The first term n\6 
is the usual decaying mode, which we are going to neglect in the following, while n 2 is 



related to the power spectrum, see e.g. p4 |. 



3. Analysis of barotropic UDM models on the pressure-density plane 

A common way to study the properties of the EoS of DE is to consider the (dw jd In a) — w 



phase space (see e.g. [74, 75, [76], f77|). Here we follow another approach, studying our 
models in the pressure-density plane, see Fig. [l]. There are several motivations for this 
choice. First of all, in the barotropic case we are considering the pressure is a function 
of the density only, so it is natural to give a graphical description on the p — p plane. 
Second, this plane gives an idea of the cosmological evolution of the dark fluid. Indeed, in 
an expanding Universe (H > 0) Eq. (E^) implies p < for a fluid satisfying the null energy 



condition [78] w > —1 during its evolution, hence there exists a one-to-one correspondence 
between time and energy density. Finally, in the adiabatic case the effective speed of 
sound we have introduced in Eq. (2^) can be written as c 2 = dp/dp, therefore it has an 



immediate geometric significance on the p — p plane as the slope of the curve describing 
the EoS p = p{p). 

For a fluid, it is quite natural to assume c 2 > 0, which then implies that the function 
p(p) is monotonic, and as such crosses the p = —p line at some point /9a- 3 From the point of 
view of the dynamics this is a crucial fact, because it implies the existence of an attracting 
fixed point (p = 0) for the conservation equation (fE3|) of our UDM fluid, i.e. p\ plays the 



2 This solution is exact if the speed of sound satisfies the equation 2c" c s — 3 (c' s ) 2 = 0, which implies 

4 

c s = 



(Cl7? + C2) 2 



where ci and C2 are generic constants. A particular case is when c\ = 0, for which the speed of sound is 
constant. 

3 Obviously, we are assuming that during the evolution the EoS allows p to become negative, actually 
violating the strong energy condition JtSJ , i.e. p < — p/3 at least for some p > 0, otherwise the fluid would 
never be able to produce an accelerated expansion. 
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Figure 1: The UDM p — p plane with the most important areas, (see the text for more detail). 
The dashed line represents the p = —p line; the dash-dotted line represents the p = —p/3 line, the 
boundary between the decelerated expansion phase of the Universe and the accelerated one; the 
dotted line p = —p/10 represents a fictitious boundary, above which the CDM-like behaviour of 
the UDM fluid dominates. The pressure and the energy density are normalised to p\. The ACDM 
model is represented here by the solid horizontal line p/pa = — 1, while the line p = represents an 
EdS model, i.e. pure CDM. 



role of an unavoidable effective cosmological constant. The Universe necessarily evolves 
toward an asymptotic de-Sitter phase, a sort of cosmic no-hair theorem (see |7^, [8(| and 
refs. therein and |6^, |54j). 

We now summarise, starting from Eqs. ( |2. 1^2.51 ) and taking also into account the 
current observational constraints and theoretical understanding, a list of the fundamental 
properties that an adiabatic UDM model has to satisfy in order to be viable. We then 
translate these properties on the p — p plane, see Fig. |l[ 

1. We assume the UDM to satisfy the weak energy condition: p ^ 0; therefore, we are 
only interested in the positive half plane. In addition, we assume that the null energy 
condition is satisfied: p + p > 0, i.e. our UDM is a standard (non-phantom) fluid. 
Finally, we assume that our UDM models admit a pa, so that an asymptotic w = — 1 
is built in. 

2. We demand a dust-like behaviour back in the past, at high energies, i.e. a negligible 
pressure p <C p for > pa- 4 In particular, for an adiabatic fluid we require that at 

4 Note that we could have p ~ —pA and yet, if p ^> pA, the Universe would still be in a matter-like era. 
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recombination \w rec \ < 10 6 , see [81, 37, 54, 55 1 



3. Let us consider a Taylor expansion of the UDM EoS p(p) about the present energy 
density po : 

p^Po + a(p-p ) , (3.1) 



i.e. an "affine" EoS model |57], 54, 55, |67|] where a is the adiabatic speed of sound 
at the present time. Clearly, these models would be represented by straight lines in 
Fig. [l], with q the slope. The ACDM model, interpreted as UDM, corresponds to the 
affine model ( |3.1| ) with a = (see [67] and [54, [37|) and thus it is represented in Fig. 
by the horizontal line p = —p\. From the matter power spectrum constraints on 
affine models [37], it turns out that a < 10~ 7 . Note therefore that, from the UDM 
perspective, today we necessarily have w ~ —0.7. 



Few comments are in order. From the points above, one could conclude that any 
adiabatic UDM model, in order to be viable, necessarily has to degenerate into the ACDM 
model, as shown in |5(| for the generalised Chaplygin gas and in [37] for the affine adiabatic 
model 5 (see [82, |83|, ^4j] for an analysis of other models). In other words, one would conclude 
that any UDM model should satisfy the condition c 2 <C 1 at all times, so that k 2 3> k 2 for 
all scales of cosmological interest, in turn giving an evolution for the gravitational potential 
$ as in Eq. 

$k * A k (l - ^ J a 2 dr^j , (3.2) 

where = <3?k (0) T m (fc), $k (0) is the primordial gravitational potential at large scales, 
set during inflation, and T m (k) is the matter transfer function, see e.g. [jS5|[. 

On the other hand, let us write down the explicit form of the Jeans wave number: 



P 



(1 + w) 



2(1 + z) 



1 



dc 2 s 3(c? 
dp 



w 



2(c s 2 



w 



6(1 +w) 



+ 



1 



(3.3) 



Clearly, we can obtain a large k 2 not only when c 2 — > 0, but also when c 2 changes rapidly, 
i.e. when the above expression is dominated by the p dc 2 /dp term. When this term is 
dominating in Eq. ( [3. 3D , we may say that the EoS is characterised by a fast transition. 

Thus, viable adiabatic UDM models can be constructed which do not require c 2 <S 1 at 
all times if the speed of sound goes through a rapid change, a fast transition period during 
which k 2 can remain large, in the sense that k 2 <C k 2 for all scales of cosmological interest 
to which the linear perturbation theory of Eq. ( [2.5[ ) applies. From point 3 above, at late 
times we must have p ~ — Pa; on the other hand, at recombination we have |tt> rec | < 10 -6 
and the speed of sound is negligible, implying p ~ constant. Therefore, the transition 
will mark the passage from a very small (possibly vanishing) almost constant p to the 



5 From the point of view of the analysis of models in the p — p plane of Fig. Ill the constraints found by 
Sandvik et al JS^] on the generalised Chaplygin gas UDM models and by Q on the affine UDM models 
simply amount to say that the curves representing these models are indistinguishable from the horizontal 
ACDM line. 
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asymptotic value p ~ ~PA or, in other words, from a pure CDM-like early phase to a post- 
transition ACDM-like late epoch. In addition, we may expect the transition to occur at 
relatively high redshifts, high enough to make the UDM model quite similar to the ACDM 
model at late times. Indeed, again from point 3 above, we infer that the fast transition 
should take place sufficiently far in the past, in particular during the dark matter epoch, 
when p ^> p\. Otherwise, we expect that it would be problematic to reproduce the current 
observations related to the UDM parameter w, for instance it would be hard to have a 
good fit of supernovae and ISW effect data. 

In the rest of the paper, in order to quantitatively investigate observational constraints 
on UDM models with fast transition, we introduce and discuss a toy model. In particular, 
we will explore which values of the parameters of this toy model fit the observed CMB and 
matter power spectra. 

Finally, let us make a last remark on building phenomenological UDM (or DE) fluid 
models intended to represent the homogeneous FLRW background and its linear pertur- 
bations. A fast transition in a fluid model could be characterised by a large value of c^, 
even larger than 1. As far as the FLRW background evolution is concerned, this fact does 
not raise any issue: the background is homogeneous, and Cg = p/p does not actually rep- 
resent a speed of sound, as there is nothing that could propagate in this case. For linear 



perturbations, at scales such that k <C kj the solution of the Eq. (2Jj) for the gravitational 
potential is Eq. (|3.2| ). Therefore, for such scales there is no superluminal propagation. 
This is because Eq. (|2.5| ) is the Fourier component of a wave equation with potential 9"/9, 
and this potential does not allow propagation for fc<ij, In building a phenomenological 
fluid model, we can therefore choose values for the parameters of the model in order to 
always satisfy the condition k kj for all k of cosmological interests to which linear theory 
applies, hence such a fluid model will be a good causal model for all scales that intends to 
represent. In other words, we can always build the fluid model in such a way that all scales 
smaller than the Jeans length A Aj correspond to those in the non-linear regime, i.e. 
scales beyond the range of applicability of the model. So, for these scales, no conclusions 
can be derived from the linear theory on the behaviour of perturbations of a UDM model 
with Cg > 1. To investigate these scales, one needs to be beyond the perturbative regime 
investigated here, possibly also increasing the sophistication of the fluid model in order to 
properly take into account the greater complexity of small scale non-linear physics and to 
maintain causality. 

4. A toy model with fast transition 

In the present section we introduce a toy model based on a hyperbolic tangent EoS, which 
is conveniently parametrised as 



P = ~PA 



tanh ( ^ 



1 - tanh 



(4.1) 



and is depicted in Fig. [2] for a typical shape. In the EoS fl4.1|) p t is the typical energy 
scale at the transition, p s is related to the rapidity of the transition, p\ is the effective 
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cosmological constant, i.e. p(p\) = — Pa- This model reduces to a ACDM, which in the 
UDM language of the previous section is described by p = — p\, in two limits: 
p t — > oo and p s — > oo. 



P 




1 2 3 4 5 6 7 



P/P A 

Figure 2: Illustrative plot of the EoS as a functions of the energy density for the hyperbolic tangent 
model. The parameters values are pt/pA = 5 and p s / Pa — 1- The energy density and the pressure 
are normalised to p\. The other five lines are the ones plotted and described in Fig. 1. 



The main properties of the EoS ( fj.l| ) are the following: 

1. The asymptotic behaviour of the pressure for p 3> pt is p ~ 0. From the considerations 
of the previous section, we expect p t 3> p\, which corresponds to a minimum value 
of the redshift z t . For instance, we have zt > 1.85 if we want to have pt > 10pA- 
In Figs. H we plot the evolution of w as a function of redshift, for pt/pA = 10 (left 
panel) and Pt/pA = 20 (right panel), for three different choices of p s / Pa- The solid line 
represents the ACDM model, while the horizontal lines respectively represent: a pure 
CDM model for w = 0; the boundary between the decelerated and the accelerated 
expansion phases of the Universe for w = —1/3. Clearly, from both panels, models 
with larger p s /pt ratio have a background evolution more similar to that of the 
ACDM model at all times. On the other hand, a smaller p s /pt ratio implies a faster 
transition between the CDMdike phase and the ACDM phase. In addition, we have 
the confirmation that the transition has to take place sufficiently far in the past, i.e. 
Pt S> pa, in order for the late time evolution of w to be in any case close to that of 
the ACDM model. 
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z z 

Figure 3: Evolution of the UDM parameter w = p/p in the hyperbolic tangent model, for pt/ PA = 
10 (left panel) and pt/pA = 20 (right panel). For reference we also plot: the w = line representing 
a flat pure CDM model (an EdS Universe); the w = — 1/3 line representing the boundary between 
the decelerated and the accelerated phases; the solid curve representing the evolution of w for a 
typical ACDM model with = 0.7. In each panel, the three dashed, dash-dotted and dotted 
lines respectively correspond to p s / Pa = E 10, 100. Clearly, the dotted lines correspond to UDM 
models with a slow transition, almost indistinguishable from a ACDM at all times, while the dashed 
lines well represent UDM models with a very fast transition from a pure CDM to a typical ACDM 
behaviour. The higher pt/pA, the earlier the UDM w transits to that of a ACDM. 



2. The speed of sound is the following: 

„ 1 tanh 2 '' 
2 PA 



Ps 1 - tanh ( 

P 



(4.2) 



illustrated in the left panel of Fig. ^. It attains its maximum value 

c!\ = ^ M..3) 

b I max / - - \ 



1 — tanh 



PA-fh 



in p = p t - For our analysis, there are two main cases to consider, assuming p t S> pa'- 

a) p t p s - In this case, | max ~ Pa/ Ps ~ 0, so that the model is close to a ACDM 
at all times. 



max 



b) pt ^> p s - In this case, we have two subcases: hi) pa < 2p s , for which Cg 

Pa/^Ps < 1 or bit) pa > 2p s , for which Cg| max ~ pA/2p s > 1. The latter subcase 
may in principle imply superluminal perturbations; fortunately, as we shall see, 
a-causal effects can be avoided if the transition is sufficiently fast. 
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1 2 345678 910 2 3 4 56789 10 

P/P A P/P A 



Figure 4: Illustrative plots of the speed of sound and p dc^/dp as functions of the energy density 
for the hyperbolic tangent model. The parameters values are pt/pA = 5 and p s /pA = 1- The energy 
density and the pressure are normalised to p\. 



As we explained in the previous section, in order to have a fast transition we must 
have pd,Cs/dp » 1 in Eq. fl3.3| ) for the Jeans wave number. This quantity is depicted 
in the right panel of Fig. |I[ For the EoS ( |4.1[ ) the derivative of is 



*£ = tonh (ez*\ 1 - famh TO = _i tanh (ez « 

d P Ps V Ps /i_tanh[^=^] P« V A 



2^f tanh ( — — ) / _ { = -— tanh ( — — ) c z s , (4.4) 

which attains its extrema at p = pt ± p s tanh -1 n/3/3) — pt ± 0.66p s . The maximum 
corresponds to the minus sign. 

Clearly, the derivative of Cg is important only in the case b) of the previous point. In 
this case the maximum is: 



^ dp 



d4 

^ dp 



. (4.5) 

Ps Ps 



For subcase bii), > 1, we always have PAPt/Ps 3> 1> while in subcase hi) there is 
also the possibility that PAPt/Ps be small. 

Let us now consider the case when the transition takes place at the lower limit 
Z\ ~ 1.85, corresponding to pt ~ 10pA- In this case, from Eq. (|4.5|), the maximum 
is Wp\/pg. Therefore, in order to have a fast transition, we must have pa > p s - 
Then, it is inevitable from point bii) that, if we want the fast transition to take place 
just before the accelerated phase of the expansion of the Universe, we must have 
> 1. In this case, shortly after the transition the pressure rapidly approaches the 
asymptotic value ~pa- 
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5. Analysis of the Jeans wave number during the transition 



The Jeans length is a crucial quantity in determining the viability of a UDM model, because 
of its effect on perturbations, which is then revealed in observables such as the CMB and 
matter power spectra. We now focus on the Jeans wave number for the toy UDM model 
introduced in the previous section and investigate its behaviour as a function of the speed 
of sound, in particular around p = p t , in the middle of the transition where the speed of 
sound is at its peak. 

Starting from the classification we presented in point 2 of the previous section, we are 
interested in the case b), namely p t 3> p s , because in this regime a fast transition in the EoS 
takes place. The majority of the adiabatic UDM models considered so far in the literature 
belongs to the case a) pt *C p B of point 2 of section f|. For the toy model Eq. (4J) as well, 
pt *C p s implies that the pressure tends to p ~ — p\ at all times, i.e. to a ACDM, as shown 
in Fig. |. 

In the case of a fast transition, from Eq. ( |3.3| ) for the Jeans wave number, it is in- 
teresting to compare the term pdc^/dp with the remaining ones contained in the squared 
brackets, namely: 

2 b{l + w) 6 

In Figs. we plot pdc^/dp, B and the Jeans wave number kj as functions of p/pa- For 
the calculation of kj we use p\ = QaPo, with £l\ = 0.7 and the critical energy density 
po = 3i?Q, where Hq is the Hubble constant. We choose pt = 100/5 a in order to consider 
a transition sufficiently back in the Dark Matter epoch (see point 1 of section and 
vary the ratio p s /pA, with p B = 10/?a, 10 _1 /3a, lCT 4 /^, in order to show examples of faster 
transitions. From the plots in Figs. it is clear that the smaller p s /pA is, the larger is 
the difference between pdc^/dp and B. 

Moreover, from Fig. [|, pdc^/dp is negative for p > p t , then for p « p t it increases 
becoming positive and intersecting the B curve a first time for p < p t . For smaller values 
of the energy density, pdc^/dp decreases again to zero, again intersecting the B curve. In 
Figs. ^-0, the same behaviour of the curves takes place and since the difference between 
the two curves is much larger, we have chosen a logarithmic scale. Therefore, the negative 
part of pdc^/dp has been omitted. 

The intersection points between the curves pdc^/dp and B represent the moments at 
which the Jeans wave number kj vanishes, as it can be seen from the right panels of Figs. 
|5||7|. In general, around these points the corresponding Jeans length becomes very large, 
possibly causing all sort of problems to perturbations, with effects on CMB and structure 
formation in the UDM model. On the other hand, for sufficiently small p s the transition 
is fast enough that i) in general the Jeans wave number becomes larger and it) it becomes 
vanishingly small for extremely short times, so that the the effects caused by its vanishing 
are sufficiently negligible, as we are going to show in the next section for the CMB and 
matter power spectra. As illustrated in Figs. HQ, by choosing progressively smaller values 
of p s we can obtain progressively larger Jeans wave numbers, while the curve kj(p) starts 
to show a plateau shape around the transition. 
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P/P A P/P A 

Figure 5: Left panel: evolution of pdc^/dp (solid line) and B (dashed line) as functions of the 
energy density. Right panel: evolution of the Jeans wave number as function of the energy density 
p/ PA- The Jeans wave number kj is in h Mpc -1 units and it has been calculated assuming fl\ = 0.7. 
The choice of the parameters is: p t — 100pA {z t ~ 5.14) and p s = 10/9A- 




Figure 6: Same as in Fig. 4, again with p t — IOO^a (zt — 5.14) but now with p s = 10 1 p\. 



Clearly, we are interested in the value of kj during the transition, because before and 
after that the negligible speed of sound implies a vanishing Jeans length, or a very large 
kj. In essence, for a fast enough transition the "average" value of kj around the transition 
is approximated by its value on the plateau - say kj - and this is, on average, the minimum 
value of kj, i.e. the maximum Jeans length for the given values of the parameters p s and 
pt- Thus, we now want to establish a relation between p s and p t for any given kj. This, 
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fixing a fcj which allows for an acceptable matter power spectrum which fits observational 
data, will help us to find the p s needed to have the transition at p t . 

The relative maximum of the Jeans wave number between the two zeros of the curve 
kj(p) corresponds approximatively to where dc^/dp assumes its maximum value, i.e. in 
p ~ p t — 0.66/? s , as we have shown in point 3 of section ^. Thus, let us define with kj 
the value of kj for p = p: as required, it is of the same order of the plateau value (see for 
example Fig. 0) of the Jeans wave number during the fast transition. 





Figure 8: The parameter p s /p\ required to obtain a given fcj (left panel) or a given E (right 
panel) vs the transition redshift z t . fcj and E are the values of the Jeans wave number and the 
efficiency at p ~ pt — 0.66p s , when dc^/dp is maximum. Left panel: fcj = 0.5, 1, 10 h Mpc -1 from 
top to bottom. Right panel: E = 10, 10 2 , 10 3 from top to bottom. The solid lines represent the 
theoretical approximations for fcj and E, the circles the numerical values, see text. 



Evaluating the analytical expression ( |3.3[ ) of fcj (p) at p under the assumption p s -C Pt 



we obtain the following approximate expression: 



k] * —B^&L 
4(l + z t ) 2 Ps 



PA / PA PA 



(5.2) 



Denning V := — 4 (1 + z t ) 2 k 2 + p% / p t and making sure that we have pdc^/dp > B for 
p = p, we can then extract from ( |5.2[) the required relation between p s and pt^ 



P + A /96^/3^+P 2 -24 

^ = ¥ £ • (5-3) 

PA 12 



In the left panel plots of Fig. g, we compare the analytical approximation (|5.3| ) with the 
numerical calculations from Eq. ( |3.3| ), for p = p and for fcj = 0.5,1,10 /i Mpc -1 , as 
functions of z t . The agreement between our analytical approximation and the numerical 
calculation is clearly very good. As we can see from the figure, if we require larger values 
of kj then p s /pA must be smaller, i.e. a faster transition is needed. On the other hand, if 
the transition takes place farther in the past, i.e. for increasing values of zt, this constraint 
is less stringent. 

Having established a good approximation for kj, we now want to determine for which 
values of pt and p s this quantity is well representative of kj around the transition, i.e. when 
we have a plateau as in Fig. |7[ In particular, this can be estimated from the difference of the 
values of pdc^/dp and B at p = p. The larger is the difference, the faster is the transition 
and the higher is the plateau effect. We therefore define the efficiency parameter E := 
p (dc^/dp) /B\ p= p which, under the assumption p t 3> p s , can be analytically approximated 
from Eq. (p|) : 

4\/3> 

E ~ - PA . (5.4) 

6(^-) +^ + 1 



v PA / PA 

From this, we obtain a new relation between p s and p t : 



Ps 1 V96y^ - 23^ - E 

p7 " 12 E ' (5 ' 5) 

In the right panel of Fig. ||, we compare the analytical approximation Q5.5| ) with the nu- 
merical calculations, for E = 10, 10 2 , 10 3 , with very good agreement. Notice that the larger 
is the efficiency E, the smaller p s /p\ must be, i.e. a faster transition is required. 6 

It is important to stress that a large efficiency is relevant in order for kj to be a more 
representative "minimum on average" value of kj during the transition, i.e. it is not a neces- 
sary condition in order to have a fast transition or a model in good agreement with observa- 
tion. This can be understood observing the multiplicative term [3/5(1 + w)} / [2(1 + z) 2 c 2 ] 
in front of the expression (|3.3| ) of the Jeans wave number kj. Indeed, for increasing values 
of pt, this term amplifies the difference [pdc^/dpj — B during the transition, giving a larger 



Assuming p s > in Eq. (|5.5|) implies E < 4\/3 {pt/ Pa)- This limit in E can be seen in the right panel 



of Fig. ^ and in Fig. |l^, where in this limit c; 



2 

8 Max 
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kj. So, one can obtain models in good agreement with observation even if the efficiency is 
low. 

Considering the plots in Figs. |||7|, where for p s = 10/9A, lO -1 /^, 10~ 4 pA we obtain 
respectively E = 1.13,597.26,692.75, in order to have a fast transition and a pronounced 
plateau, we infer that E > 700 is needed. In addition, this requirement was also confirmed 
by the study of the matter power spectrum, see the next section. 




500 1000 



Figure 9: Efficiency E as function of z t , for fej = 5, 2, 1, 0.5 h Mpc 1 (from top to bottom). 



Substituting Eq. ( |5.3|) in Eq. (5.4), we can now obtain an approximate expression of 
£ as a function of the redshift of the transition zt, which allows us to estimate the range 
of zt (and thus pt) for which the efficiency E is above a certain threshold, for a given kj. 
In Fig. |9] we show the plot E vs z t for kj = 5, 2, 1, 0.5 h Mpc -1 (from top to bottom). For 
increasing values of kj, the range of zt in which E > 700 becomes larger, as expected. 

We can now use the relation ( |5.3D to understand how large the speed of sound can be 
during the transition. To this purpose, we substitute Eq. (5.3) in the maximum value of 



Cg, at p = pt- We plot this c 2 Max in Fig. ^ as a function of z t for kj = 10, 1, 0.5 h Mpc -1 
(from top to bottom). On the same figure we also plot the curve of the maximum of the 
speed of sound fixing the value of the efficiency at E = 700 (solid red line). Then, the 
vertical dashed line corresponds to the value of zt for which p s — > in Eq. ( |5.5| ), giving 
r 2 

SMax 



oo in Eq. (O). In order to have E > 700, we must consider the area above the 



solid red line. We can see that, for increasing values of z t , the value of Cg Max required to 
have a fixed kj decreases. For example, for kj = 1 h Mpc -1 , in order to have Cg Max < 1, 



the transition has to take place at z% > 50, while if z t ~ 5 we have Cg Max ~ 



10 s 



6. The CMB and matter Power spectra: toy model predictions 

In order to compare the predictions of our toy UDM model with observational data, we 
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C sMax 102 




500 1000 



Figure 10: Evolution of the maximum value of Cg as function of the redshift z t . The choice of p s 
as function of pt is given by Eq. ( |5.3[ ) in the text for fcj = 10, 1,0.5 h Mpc -1 (from top to bottom, 
dashed, dash-dotted and dotted line) and by Eq. (5.4) for E = 700 (solid red line). Then, the 
vertical dashed line corresponds to the value of zt for which p s — > in Eq. ( |5.5|) , giving Cg Max — > oo 



in Eq. (4.3) 



have used a properly modified version of CAMB 7 J66| for the computation of the CMB 
and the matter power spectra. In particular, we have modified the original definition of 
the density contrast for the case of adiabatic UDM models. Indeed, we have to define the 
UDM density contrast as 5 := 5p/p m [37], where here p m = p — p\ is the "aether" part of 

In this case, starting from the perturbation theory we outlined in 



the UDM fluid |67], | 

section |2[ we can infer the link between the density contrast and the gravitational potential 
via the Poisson equation in the following way: 



5 (k; z) 



5p(k;z) 



3tf 2 OnO (1 + z) 



-k : 



(3/2)0 



m0 



(6.1) 



for scales smaller than the cosmological horizon and z < z rec , where z mc is the recombina- 
tion redshift (z Tec ~ 10 3 ). 

We compare the theoretical predictions of our toy model with the WMAP 5-year 
data [18, 20] and the luminous red galaxies power spectrum measured by the SDSS 
collaboration [fT5| . The CMB data used in our plots are available on the LAMBDA 8 
website, while those regarding the matter power spectrum are implemented in a modified 
version of the CosmoMC software 9 . 



7 http:/ /camb.info/ 

8 http:/ /lambda.gsfc. nasa.gov 

9 http : // cosmologist .info/ cosmomc / 
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We consider as reference the flat ACDM model described by the best-fit parameters 
found by combining WMAP5 data with measurements of Type la supernovae and Baryon 
Acoustic Oscillations [18, 20|, with values provided on the LAMBDA website (68% CL 
uncertainties): n h0 h 2 = 0.02265 ± 0.00059, n m0 h 2 = 0.1143 ± 0.0034, 0\ = 0.721 ± 0.015, 
H = 70.1 ± 1.3 km s^ 1 Mpc" 1 , n s = 0.960±ft|g£ and A| = (2.457±g:g^) x 10" 9 . For our 
toy model, we keep the same amount of baryons but choose a vanishing CDM content. 

In Figs. |TT|-[r4l we plot the theoretical predictions of our model, those of the reference 
ACDM and the observed CMB and matter power spectra data. Each of Figs. 11-14] re- 
spectively correspond to pt/pA = 10 5 , 10 3 , 10 2 , 10, i.e. to zt ^ 66, 13, 5.7, 2. Guided by the 
analysis in section [51 for each transition density pt we have chosen values of p s which clearly 
show the progressive enhancement of the agreement between the predicted matter power 
spectrum and the observed one. Moreover, in the matter power spectrum plots, for each 
choice (pt)Ps) we draw a vertical dashed line representing the corresponding value of kj. 

We can see from Figs. 11 and 12 that the CMB anisotropies predicted by the reference 
ACDM and by our toy model are indistinguishable for a large range of p s . However, while 
at the higher transition redshift of Fig. 11 the matter power spectrum also allows the same 
broad range of p s values, at the smaller z± of Fig. 12 we start to see the need for a faster 
transition, i.e. a smaller p s , to have an acceptable power spectrum. As expected from the 
analysis of section || the effect becomes more pronounced as the transition occurs at the 
smaller and smaller redshifts of Figs. 13 and 14. In the case p s /PA ~ 0.1 of Fig. H the first 
acoustic peak of CMB is higher with respect to the observational data. This effect can be 
explained by looking at the matter power spectrum. Indeed, the latter moves away from 
the reference ACDM before the equivalence wavenumber k eq ~ 0.01 h Mpc -1 . In other 
words, the gravitational potential starts to oscillate and to decay for k < k cq , therefore 
affecting those modes entering the horizon before the matter-radiation equivalence epoch. 

Finally, in Fig. [14| the first CMB spectrum peak is lower than the observed one for any 
value of p s . Note in the left panel of Fig. [3| that, for p t = WpA, the background evolution 
is sensibly different from the reference ACDM. Indeed, in this case our model behaves 
like a pure CDM Einstein-de Sitter for a much longer time. One possibility to avoid this 
discrepancy is to slightly increase Qa- Therefore, again using p t = 10pA and p s = 10~ 5 pA, 
in Fig. 15 we have chosen ^a = 0.721, 0.742, 0.772, with the first value again corresponding 
to the above mentioned reference ACDM and the other two corresponding to the best-fit 
and its upper uncertainty obtained by WMAP5 using CMB data only (see Jl9| and the 
LAMBDA website). The agreement between the CMB prediction and the observational 
data is again good for ^a = 0.742, with a good matter power spectrum. 



7. Conclusions 



The last decade of observations of large scale structure [14, [H| 16, 17, 21, the search 



for la supernovae (SNIa) 0, |8[ ||, [l(J and the measurements of the CMB anisotropies 
11, [l8|, |l9], are very well explained by assuming that two dark components govern the 
dynamics of the Universe. They are DM, thought to be the main responsible for structure 
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formation, and an additional DE component that is supposed to drive the measured cosmic 
acceleration 11, 12, p^| . 

At the same time, in the context of General Relativity, it is very interesting to study 
possible alternatives. A popular one is that of an interaction between DM and DE, without 
violating current observational constraints 11, [28|, EHJ, |30|, 31, p2l B3|. This possibility could 



alleviate the so called "coincidence problem" |23], namely, why are the energy densities of 
the two dark components of the same order of magnitude today. Another attractive, 
albeit radical, explanation of the observed cosmic acceleration and structure formation is 



to assume the existence of a single dark component: UDM models |38|, |39|, |50|, |48|, |49|, ^7 



52| , p8| , |3q , p3| , |57| , [55| , |62| , p4| , p7| , |65(1 where, by definition, there is no coincidence problem. 



In the present paper we have investigated the general properties of UDM fluid models 
where the pressure and the energy density are linked by a barotropic equation of state 
(EoS) p = p(p) and the perturbations are adiabatic. Using the pressure-density plane, 
we have analysed the properties that a general barotropic UDM model has to fulfil in 
order to be viable. We have assumed that the EoS of UDM models admits a future 
attractor which acts as an effective cosmological constant, while asymptotically in the past 
the pressure is negligible, studying the possibility of constructing adiabatic UDM models 
where the Jeans length is very small, even when the speed of sound c s is not negligible. 
In particular, we have focused on models that admit an effective cosmological constant 
and that are characterised by a short period during which the effective speed of sound 
varies significantly from zero. This allows a fast transition between an early epoch that is 
indistinguishable from a standard matter dominated era, i.e. an Einstein de Sitter model, 
and a more recent epoch whose dynamics, background and perturbative, are very close to 
that of a standard ACDM model. 

In the second part of the paper, in order to quantitatively investigate observational 
constraints on UDM models with fast transition, we have introduced and discussed a toy 
model based on a hyperbolic tangent EoS [see Eq. ([4.1])]. We have shown that if the 
transition takes place early enough, at a redshift zt > 2 when the effective cosmological 
constant is still subdominant, being also fast enough, then these models can avoid the 
oscillating and decaying time evolution of the gravitational potential that in many UDM 
models causes CMB and matter fluctuations incompatible with observations. Consequently, 
the background evolution, the CMB anisotropy and the linear matter power spectrum 
predicted by our model do not display significant differences from those computed in a 
reference ACDM (l^, 20], because the Jeans length Aj = a/fcj, where /cj is the Jeans wave 
number [see Eq. ( |3.3[ )1, remains small at all times, except for negligibly short periods, even 
if during the fast transition the speed of sound can be large. In this way, our toy models 
(and more in general UDM models with a similar fast transition) can evade the "no-go 
theorem" of Sandvik et al ]50| , as we discussed in the introduction. 

Specifically, we have analysed the properties of perturbations in our toy model, focusing 
on the the evolution of the effective speed of sound and that of the Jeans length during the 
transition. In this way, we have been able to set theoretical constraints on the parameters 
of the model, predicting sufficient conditions for the model to be viable. Finally, guided 
by these predictions and using the CAMB code [36|, we have computed the CMB and the 
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matter power spectra showing that our toy model, for a wide range of parameters values, 
fits observation. 

The full likelihood analysis for this model and its parameters would be an interesting 
extension of the study carried out here, which we will address in a future work. 
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Figure 11: Matter power spectrum (upper panel) and CMB power spectrum (lower panel) for 
p t = 10 5 pa, i.e. a transition at z t ~ 66. The values of p s /pA are: p s /pA — 0.25 (dashed green line), 
Ps/pa — 1 (dash-dotted red line) and p s / Pa — 10 (solid blue line). The choice of the parameter 
has been done in order to have: fcj ~ 2 h Mpc -1 (dashed green vertical line) and E ~ 4 10 5 for 
Ps/pa — 0.25; fcj ~ 1 h Mpc -1 (dashed red vertical line) and E ~ 9 10 4 for p s / Pa — 1; fcj ~ 0.3 h 
Mpc -1 (dashed blue vertical line) and E ~ 10 3 for p s / Pa — 10. The reference ACDM (see text) is 
represented by the solid black line. Notice that the theoretical curves representing the CMB power 
spectrum for our models and the reference ACDM are indistinguishable. 
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Figure 12: Matter power spectrum (upper panel) and CMB power spectrum (lower panel) for 
Pt = IO'Va, i.e. z t ~ 13. The values of p s /pA are: p s / Pa — 10~ 3 (dashed green line), p s /pA — 0.05 
(dash-dotted red line) and p s /pA — 0.5 (solid blue line). The choice of the parameter has been done 
in order to have: k,j ~ 1.5 h Mpc -1 (dashed green vertical line) and E ~ 7 10 3 for p s /pA — 10~ 3 ; 
kj ~ 0.2 h Mpc" 1 (dashed red vertical line) and E - 7 10 3 for p s /p A ^ 0.05; kj - 0.07 h 
Mpc" 1 (dashed blue vertical line) and E ~ 2 10 3 for p s / Pa — 0.5. Also in this case z% ~ 13 the 
theoretical curves representing the CMB power spectrum for our models and the reference ACDM 
are indistinguishable. However, the matter power spectrum requires a faster transition, i.e. smaller 
p s values. 
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Figure 13: Matter power spectrum (upper panel) and CMB power spectrum (lower panel) for 
pt = 10 2 pa, i.e. z t ~ 5.7. The values of p s /pA are: p s /pA — 10~ 5 (dashed green line), p s /p\ — 0.01 
(dash-dotted red line) and p s /pA — 0.1 (solid blue line). The choice of the parameter has been done 
in order to have: kj ~ 3.3 h Mpc -1 (dashed green vertical line) and E ~ 6.9 10 2 for p s / p\ ~ 10~ 5 ; 
k 3 - 0.1 h Mpc- 1 (dashed red vertical line) and E - 6.9 10 2 for p s /p A ~ 0.01; fcj - 0.03 h Mpc" 1 
(dashed blue vertical line) and E ~ 6 10 2 for p s /pA — 0.1. At the relatively small transition redshift 
z t ~ 5.7 a viable matter power spectrum requires an even faster transition, i.e. smaller p s values. 
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Figure 14: Matter power spectrum (upper panel) and CMB power spectrum (lower panel) for 
pt = 10pA, i-e. zt ~ 2. The values of p s / pA are: p s /pA — 10 -7 (dashed green line), p s /p\ — 10~ 5 
(dash-dotted red line) and p s /p\ — 10~ 4 (solid blue line). The choice of the parameter has been 
done in order to have: fcj ~ 7 h Mpc -1 (dashed green vertical line) and E ~ 70 for p s / Pa — 10~ 7 ; 
fej - 0.7 h Mpc -1 (dashed red vertical line) and E ~ 68 for p s /p A ~ 10~ 5 ; fcj - 0.23 /i Mpc -1 
(dashed blue vertical line) and E ~ 68 for p s / pa — 10 -4 . At this smaller transition rcdshift 
z t ~ 2 the background evolution in our model is so strongly modified that, no matter how fast the 
transition is, i.e. even for p s values giving an acceptable matter power spectrum, there is no way 
to fit the CMB first peak for the given choice of ft\, i.e. the same than the reference ACDM (see 
text). 
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Figure 15: Matter power spectrum (upper panel) and CMB power spectrum (lower panel) for one 
of the transition model of Fig. |l4|, i.e. pt = Wp\ and p s /pA — 10~ 7 , with a transition at Zt ~ 2, this 
time for different values of 51a- The solid blue, red dot-dashed and green dashed lines respectively 
correspond to fi A = 0.721,0.742,0.772. The reference ACDM (solid black line) is the same of the 
other figures. With a slightly higher ft A = 0.742 (WMAP5 best fit value with CMB data alone 
|l9||) our model now produces an acceptable fit. 
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